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Holographic melting and related properties of mesons in a quark gluon plasma 
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We analyse mesons at finite temperature in a chiral, confining string dual. The temperature 
dependence of low-spin as well as high-spin meson masses is shown to exhibit a pattern familiar 
from the lattice. Furthermore, we find the dissociation temperature of mesons as a function of their 
spin, showing that at a fixed quark mass, mesons with larger spins dissociate at lower temperatures. 
The Goldstone bosons associated with chiral symmetry breaking are shown to disappear above the 
chiral symmetry restoration temperature. Finally, we show that holographic consideration imply 
that large-spin mesons do not experience drag effects when moving through the quark gluon plasma. 
They do, however, have a maximum velocity for fixed spin, beyond which they dissociate. 
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I. INTRODUCTION AND SUMMARY 

Of the various gravity duals to confining, non- 
supersymmetric gauge theories in four dimensions, the 
Sakai-Sugimoto model [l|, [|J stands out because it con- 
tains chiral fermions as well as an associated chiral sym- 
metry breaking transition. This model has been anal- 
ysed in quite some detail at zero temperature. More re- 
cent work has analysed the phase diagram of the model, 
showing an interesting structure in which the confine- 
ment/deconfinement transition does not necessarily co- 
incide with the chiral symmetry breaking transition Q 
(see also In the present paper, we will analyse the 
various phases which were found in Q in more detail, 
focussing on the properties of mesons. 

The D4-brane background of the Sakai-Sugimoto 
model with Euclidean signature exhibits two circle direc- 
tions. The analysis of [3( has shown that this geometrical 
structure, plus the presence of the D8-probes, leads to 
three different phases (when the ratio L/R, see figure [TJ 
is sufficiently small). In the low-temperature phase the 
background is the Euclidean continuation of a Lorentzian 
background without any horizon. Gluons are confined in 
this phase. After the confinement/deconfincment transi- 
tion for the gluons, we enter the intermediate-temperature 
phase, in which the Euclidean geometry is instead given 
by the continuation of a black brane geometry. In this 
phase gluons are deconfined, but chiral symmetry is still 
broken. Mesonic bound states still exist, as the D8-brane 
embedding is not yet touching the horizon. At sufficiently 
high temperature, the lowest-energy configuration of the 
D8-branes is one in which they hang vertically from in- 
finity down to the horizon. This is the high-temperature 



* Email: kasper. peeters@aei.mpg.de 
t Email marija.zamaklar@aei.mpg.de 
t Email cobi@pos t.tau.ac.il| 



phase, in which chiral symmetry is restored. If the ra- 
tio L/R > 0.97, there is no intermediate-temperature 
phase, so that the confinement/deconfincment and the 
chiral symmetry breaking transition coincide. For more 
details, we refer to 

In the present paper, we will have a close look at 
the spectrum of low-spin as well as high-spin mesons in 
these three phases. From both lattice results and exper- 
iment, it is known that the meson spectrum shows inter- 
esting behaviour as the temperature goes up. There is 
abundant evidence, for instance, that heavy-quark bound 
states can survive in a quark-gluon plasma to tempera- 
tures which are substantially higher than the confinc- 
ment/deconfinement transition [5[. Moreover, heavy- 
quark mesons tend to survive to higher temperatures 
than light-quark mesons (see e.g. @, 0])- We would like 
to see to which extent the Sakai-Sugimoto model is able 
to say something qualitatively about such situations. 

Without further ado, let us briefly summarise the out- 
come of our analysis. Firstly, we find that the tempera- 
ture dependence of both small-spin (figure 2]) and large- 
spin (figure |T3J) mesons is qualitatively similar to what is 
known from mesons in a quark-gluon plasma pi. \@L [ic| . 
That is, the masses of mesons with low excitation num- 
ber go down as the temperature goes up. The p and ai 
meson masses approach each other, although they do not 
become degenerate at the temperature where chiral sym- 
metry is restored. For higher excitation level, the temper- 
ature dependence appears to be more complicated (fig- 
ure o. 

Secondly, we analyse the temperature dependence of 
large-spin mesons, which are modelled by rotating U- 
shaped strings hanging from the probe brane. We find 
that for large-spin mesons there is a maximum value of 
angular momentum beyond which mesons cannot exist 
and have to "melt", i.e. dissociate (figure fTTj) . For a 
given spin, the maximum temperature can be substan- 
tially lower than the temperature of the chiral symmetry 
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FIG. 1: The Sakai & Sugimoto model in the low-temperature (left panel), intermediate-temperature (middle panel) and high- 
temperature (right panel) phases. Depicted are the two compact Euclidean directions x 4 (the Kaluza-Klein circle) and t (Wick 
rotated time). The D8 embedding is shown as well (in blue). 



restoration transition. We consider this maximum-spin 
effect a prediction of the string/gauge-based models, as 
it is quite generic and qualitatively independent of the 
specifics of the background. 

Thirdly, we find that large-spin mesons do not expe- 
rience a drag force. Instead, there is a critical velocity 
beyond which a fixed-spin meson will dissociate. Equiv- 
alently, the maximum interquark distance goes down as 
the velocity is increased (figure [T6|) . The absence of drag 
is unlike the situation for qua rk p robes in a quark-gluon 
plasma QH GJ, [H [3 G3, M, Hf • For large-spin mesons 
at finite temperature, one can find generalised solutions 
where the meson moves with linear velocity, rigidly, with 
free boundary conditions in the direction of motion (fig- 
ure IT5|) . Hence one does not need to apply any force 
to maintain this motion. In the dual language, this is a 
reflection of the fact that if the quark gluon plasma is 
not hot enough to dissociate mesons, then these colour 
singlets will not experience a drag force generated by 
monopole interactions with the medium. 1 



trivial metric, 



II. REVIEW OF THE SYSTEM AT VARIOUS 
TEMPERATURES 



In this section we review the various phases of the 
Sakai-Sugimoto model, originally proposed in [l|, Q and 
analysed at finite-temperature in Q. Let us first recall 
the Euclideanised model at zero temperature. The pure 
glue theory is described by a background with a non- 



1 The interaction with the medium could be still via dipole (or 
higher moment) interactions. These, however, require coupling 
of the string sigma model to the fields on the worldvolume of the 
brane, and will be discussed elsewhere. 
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as well as a non-constant dilaton and a four-form RR 
field strength given by 

3/4 



e+ = g 



n D4 



i?D4 

7rg s NJ 3 s . 



2ttN c 

F4 = ^T £4 



(2) 



Here u is the radial direction, which is bounded from be- 
low by u > ua- Our four-dimensional world is along the t 
and x % (i = 1,2,3) directions. The main characteristic 
feature of the background is the submanifold spanned 
by x 4 and u, which has the geometry of a cigar. The tip 
of the cigar is non-singular if and only if the periodicity 
of x 4 is 
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Therefore, the circle x 4 shrinks from its asymptotic 
size 27ri? down to zero and smoothly caps off at u = u A - 
Chiral quarks in the fundamental representation are in- 
corporated by introducing a probe D8-brane into the sys- 
tem. The probe fills out the full space, except in the sub- 
manifold of the cigar, where it has a non-trivial profile. 
At zero temperature, the spectrum of the glueballs (ob- 
tained from the fluctuations of the supergravity fields of 
the background), as well as of low-spin mesons (obtained 
from the fluctuations of the fields on the probe brane) is 
discrete, signalling that we are in the confining phase. 

Turning on the temperature T corresponds to the com- 
pactification of the Euclidean time direction t. The tem- 
perature is related to the size of the compact time direc- 
tion through the standard field theory relation (3 = 1/T. 
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We now have two circle directions, one being the asymp- 
totic circle of the cigar discussed above, with period 2ttR, 
and the other one being the thermal circle, with period (3. 
This geometry is depicted in the left panel of figure [TJ 
which also shows the embedding of the D8-branc. The 
glueball and meson spectra in this phase are still dis- 
crete. In fact, they are unmodified with respect to the 
zero temperature situation, given that the metric is un- 
modified (modulo the global identification) . We will refer 
to this phase as the low-temperature phase. 

As the temperature of the system is increased, the sys- 
tem will at a certain point undergo a first order phase 
transition in which the gluonic degrees of freedom get 
deconfined. This is reflected in the change of the back- 
ground metric from (TTJ) to 
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In the two backgrounds (JTJ) and ([4} , the DBI actions have 
the following form, 
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where from now on we will drop the subscript on / as it 
is clear from the context which one should be used. In 
the confining background, the solution to the equations 
of motion (depicted by the blue curve on the left panel 
of figure [1]) is given by 
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The metric has the same form as the low-temperature 
metric ([1]), but with the role of the t and x directions ex- 
changed. This is depicted in the middle and right panels 
of figure [TJ The phase transition between the metric (TTJ) 
and the metric (U happens when their free energies are 
equal, which is when the two circles have the same radius, 
i.e. when T = T c = l/2nR. This transition is first-order 
because the solutions do not smoothly connect, and con- 
tinue to exist as separate solutions both below and above 
the transition. The spectrum of the glueball fluctuations 
in the phase described by (jU) is continuous, signalling the 
deconfinement of the gluonic degrees of freedom. 

An order parameter of this phase transition is the 
Polyakov loop, corresponding to a string wrapped around 
the time direction. In the low-temperature phase the 
time circle is non-contractible and not the boundary of 
a disc, so that the Polyakov loop vanishes. After the 
transition, this circle becomes contractible, resulting in a 
non-zero expectation value for the Polyakov loop. Other 
order parameters are the Wilson loop (which has a linear 
quark/anti-quark potential in the first background but 
vanishing tension in the second) and the behaviour of 
the free energy as a function of N c (namely ~ N® for the 
first background and ~ N 2 for the second). 

In contrast to the gluonic degrees of freedom, the mat- 
ter degrees of freedom may still have bound states in this 
phase. In order to analyse the matter phases, one needs 
to determine the shape of the D8-brane, which can be 
obtained from the DBI part of the effective action, 
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The D8-brane solution in the middle panel is related to 
this one by the simple expression 
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However, there is a second configuration in the de- 
confined background, in which the embedding is simply 
given by x A = const, (depicted in the third panel of fig- 
ure!]]). The most important fact about this third phase is 
that chiral symmetry is restored. Because there are two 
stacks of branes, every mode will appear twice, as a rep- 
resentation of lJ(Nf)i J and of XJ(Nf)n. So we potentially 
have two phases in the deconfined background, which we 
will call the intermediate- and /lig/i-temperature phase. 
Whether or not the intermediate-temperature phase is 
realised depends on the ratio L/R, where L is the asymp- 
totic distance between the D8 and anti-D8 (see fig- 
ure [JJ). In H it was found that for L/R > 0.97, the 
intermediate phase is absent altogether, i.e. the con- 
finement /deconfinement and chiral symmetry restoration 
transitions occur simultaneously. One can map the pa- 
rameter L/R to a gauge dynamical parameter by noting 
that 1/R is the scale of the glueball masses, while (as we 
will show in section ITTT)) . l/L is proportional to the mass 
of mesons composed of quarks of large constituent mass, 
ito ^> ut- We refer the reader to Q for more details on 
the full phase diagram, including the region where the 
supergravity picture is not valid. 

In the string/gauge duality models, one associates 
low-spin mesons to fluctuations of the gauge fields and 
(pseudo) scalars that live on the probe branes, while high- 
spin ones are associated to spinning string configurations. 
We will address the thermal spectrum of the former in 
section ITTTl and the latter in llVI 
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III. LOW-SPIN MESONS AT INTERMEDIATE 
TEMPERATURE 

In this section we will analyse the spectrum of low-spin 
mesons (with spin < 1) in the intermediate-temperature 
regime. The spectrum of low spin mesons in the low- 
temperature phase is unmodified with respect to zero 
temperature. This is a consequence of the fact that the 
(Euclidean) metric is globally unmodified (the only differ- 
ence is that there is now a global periodic identification of 
the Euclidean time direction). This perhaps unexpected 
feature seems to be a generic property of large- N c theo- 
ries SGilGJ] (see also section Q1TX]) . 

In the intermediate-temperature phase we expect that 
the spectrum of low-spin mesons is discrete, because the 
probe does not intersect the horizon. This is similar to 
the mechanism which ensures discreteness of the spec- 
trum at zero temperature. Furthermore, given that the 
effective tension of strings near the brane decreases with 
the increase of temperature, one expects that the masses 
of the mesons decrease as the probe brane comes closer 
to the horizon, or equivalently, as the temperature is in- 
creased. Indeed, our explicit computation shows this be- 
haviour. 

Our starting point is the metric of the background de- 
scribing the hot gluonic plasma, as given in (g]). This 
background leads to an induced metric on the D8-brane 
world- volume, which reads 
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We are interested in computing the spectrum of vector 
mesons, by considering small fluctuations on the world- 
volume gauge fields of the probe D8-brane (scalar and 
pseudoscalar mesons can be treated with similar meth- 
ods but will not be discussed here). In order to do this, 
we first expand the gauge field on the worldvolume of the 
D8-brane as [l[ 
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F^u = ]T <t> n {u) - B^dMu) 
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where the last line is obtained by taking (f>t n ) = 
m~ 1 d u '4'(n)( u )- To simplify the consideration, we choose 
to focus on the space-like components Bi of the vector 
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FIG. 2: The asymptotic distance L between the D8 and anti- 
D8 as a function of tjt ■= ut/uq. The red dashed line in- 
dicates the phase transition to the high-temperature phase, 
above which one should instead use two parallel D8 stacks. 
Note that "L — const. — > uq — const." holds to good approx- 
imation in the intermediate-temperature regime. 



fields. Moreover, we only consider the masses defined by 

(11) 



dlB^ = -mlB^ 



that is, we consider the behaviour of the pole mass rather 
than the screening mass. Thus, we consider only spatially 
homogeneous modes, i.e. we consider the equation of mo- 
tion for fields satisfying 9,5^™ = 0. 2 In this case, after 
using the equations explicit probe brane embedding ([5]) 
and ((?]) the probe brane action reduces to 
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After a partial integration with respect to the u- 



(10) coordinate, the equation of motion for the field B 



(m) 



2 For non-homogeneous modes the definition of mass becomes 
more subtle. However, since we expect that the homogeneous 
modes considered here are rich enough to demonstrate all rele- 
vant features of the spectrum as a function of the temperature, 
we restrict our attention to this case. 
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FIG. 3: Meson masses (squared) as a function of level n for the zero temperature (upper dots) and intermediate-temperature 
phases (lower dots), for an asymptotic separation between the D8 and anti-D8 fixed at L/R = 0.63, with ua — ut = 1. 
(Extensions of these curves to even higher level show m 2 
which is one of the unrealistic features of these models [20l. l2ll. [2^|V 
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This equation will reduce to the canonical form (with 
thermal masses m n as in (jlip ) if the modes ipi n ) satisfy 
the equation 

- «V2 7 -V2 /(tt) V2 du ( u 5/ V 1/2 /M 1/2 a^(„)) 

= -Rd4 m 2 n V>(n) ■ (14) 

Equation (|14|) is very similar to the equation in the zero 
temperature case. The only difference is the appearance 
of extra factors f(u) x / 2 in the term on the left-hand side. 
In order to have canonically normalised kinetic terms, the 
modes should also satisfy the normalisation conditions 

duU 4 7 1/2 / _1/2 (w) ( J 1p( m )1p(n) = $mn , 



3 / p \ 3 / 2 
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The zero mode 4>^ — u 5 ^ 2 f( u ) 1 / 2 7 1 ^ 2 is normalisable 
with this norm (there is no problem at the horizon be- 
cause wo > ut), and so we see that there is still a massless 
pion present in the intermediate-temperature phase. 

Before applying numerical methods to compute the 
spectrum, it is interesting to note that in the limit of 
i*o S> ut the spectrum simplifies and one can easily de- 
termine the scale of the meson masses. In this limit, 



which corresponds to a small separation distance between 
the stacks of branes and anti-branes L R, the thermal 
factor f(u) — ► 1 and in particular also /(urj) — * 1- There- 
fore, 
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where the dimensionless quantity y = u/uq. In fact, for 
this case we can rewrite (|14[) in terms of y in the following 
form 



y /2 7 
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Now since the left-hand side is expressed in terms of the 
dimensionless quantity y, the right-hand side should also 
be dimensionless which implies that 
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Now since uq ~ 1/L 2 (see below), the final conclusion is 
that the mass of these "short" mesons scales as 



1 



(19) 



To find the remainder of the spectrum, we can solve 
equation (|14p numerically using a shooting technique. 
We will keep the asymptotic separation L between the 
D8 and anti-D8 stacks fixed, and use the expression 



L= dx 4 = 



' du 



(20) 



to determine uq for a given temperature (see also fig- 
ure O. The results for the thermal masses m 2 of the 
vector mesons are depicted in figure [3] The spectrum is 
discrete, labelled by the mode number n and parity. By 
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FIG. 4: Squared masses of the two lightest vector mesons 
as a function of temperature in the intermediate-temperature 
regime, for L = 0.256 < 0.97 R. Masses are normalised as 
in (|14p . The decrease of the mass is in qualitative agreement 
with results from other models and experiments 0], though 
one should be cautious since the latter are not in the large- 7V C 
limit. The p and ai do not become degenerate in our model. 



FIG. 5: Behaviour of six lowest-mass vector mesons (for 
fixed L/R) as a function of temperature in the intermediate- 
temperature regime, and a comparison with the masses in the 
low-temperature phase. Note the jump in the spectrum at the 
confinement/deconfinement transition, and also note that the 
masses of the higher excited modes do not exhibit monotonic 
behaviour as a function of temperature. 



comparing with the zero-temperature result, we observe 
that the masses of light mesons decrease as the temper- 
ature is increased. The temperature dependence of the 
masses of the "p" and "ai" mesons are shown in figure 2] 

Our findings are in qualitative agreement with lattice 
computations [1, [23[ and the results of the "phenomeno- 
logical" model of Karch et al. [24]], although we do not 
see the degeneracy of the p and a\ meson masses at high 
temperature. Moreover, we will see in section [V] that the 
spectrum becomes continuous in the phase where chiral 
symmetry is restored. Therefore, the phase transition 
from intermediate to high temperature, which is first- 
order in our model, is rather different from QCD. How- 
ever, what we see is consistent with the fact that mesons 
in the chirally symmetric phase are extremely unstable in 
the large- N c limit: their decay widths scale as a positive 
power of iV c . 3 

This qualitative behaviour of the spectrum in our 
model (and the ad-hoc model of [Hj]) is a direct conse- 
quence of the fact that the flavour branes are, for fixed L, 
closer to the horizon when the temperature is higher. 
Hence, this decrease of masses with temperature will be 
a common feature for all gravitational backgrounds which 
contain a horizon. However, the precise behaviour of the 
masses as a function of temperature is clearly model de- 
pendent. In our case, it depends on parity, the level n 
and the spin of the meson, as well as on the constituent 
quark masses. It is interesting to note that generically, 
for fixed spin, level and parity, the meson mass runs 
slower for larger constituent quark masses than for lighter 
quarks. This behaviour is again in qualitative agreement 



with lattice simulations, see for example [26|, [27], HH, Hil 
for quenched results and [3(| for a confirmation in un- 
quenched QCD. There, it is found that deeply bound 
heavy quark states, such as the r\ c (J — 0) or the J/ty 
( J = 1), survive up to temperatures which are < 2T cr j t , 
where T cr i t is the temperature where the light mesons 
p, u> and (f> melt. For our system, this behaviour is a nat- 
ural consequence of the fact that the constituent quark 
mass is, roughly speaking, related to the distance of the 
tip of the probe brane to the horizon. If this distance is 
increased (while keeping the horizon temperature fixed 
in units of the radius of the circle x 4 ), a meson of the 
same spin, parity and level will correspond to an excita- 
tion of the brane which is further away from the horizon, 
and hence less affected by the temperature. Note that 
because the phenomenological model of [24[ lacks probe 
branes, it does not have this parameter at its disposal. 
Hence, it thus seems harder to reproduce this feature of 
the mesonic spectrum in such ad-hoc models. 

One can also follow the mass trajectories for the higher 
excited modes. This results in figure [5] This plot shows 
interesting qualitative similarities with results in real- 
world QCD, see e.g. [9(, although the latter results are 
not in the large- iV c limit so one should be careful with 
a direct comparison. It would be very interesting to ex- 
tract more quantitative dependence of the meson masses 
on the aforementioned set of parameters from our model. 
We leave this question for future investigation. 



A. Phase transition from low to intermediate 
temperature 



For results of other, non-holographic models see e.g. the summary 
in [il. 



As mentioned before, the mesonic spectrum in the low- 
temperature phase is unchanged as the temperature is 
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increased. Although this may sound counterintuitive, it 
seems to be a generic feature of large- N c theories [H, . 
Namely, it was shown in [HI, [III] that, for example, the 
value of the chiral condensate (tptp) is independent of the 
temperature, essentially because in the confining phase 
the theory behaves effectively as a gas of non-interacting 
glueballs and mesons. This behaviour is indeed what is 
reproduced from a generic supergravity solution, since 
the metric for finite, low temperatures is the same as the 
metric for zero temperature (note that the same remains 
true even if Nf ~ N c and one constructs the back-reacted 
geometry). To incorporate finite-temperature effects on 
the meson spectrum and the value of the chiral con- 
densate, one would need to go beyond the supergravity 
approximation and deal with genuine finite-temperature 
string theory. 

Let us nevertheless see how the intermediate- 
temperature and low-temperature phases are connected 
in our model. To do this, we send the temperature to 
the critical value T c where the gluonic degrees of free- 
dom get deconfined (i.e. T — > 1/2ttR). Figure [5] shows 
that the mesonic spectrum at intermediate temperature 
is not connected to the spectrum in the low-temperature 
regime. However, since as argued in [3] the transition 
is first-order, such a jump should be expected. Do no- 
tice, however, that the mesons are only mildly affected 
by the "violent" transition of the gluonic degrees of free- 
dom. The mesonic spectrum, though discontinuous, is 
not rearranged: for every meson state before the transi- 
tion, there is a corresponding state after the transition. 
In a certain sense the mesons are merely spectators in 
this phase transition. 



IV. HIGH-SPIN MESONS AT INTERMEDIATE 
TEMPERATURE 

In this section we examine the properties of high-spin 
mesons as a function of the temperature. For low-spin 
mesons, we have seen that the dependence of the string 
tension on the temperature determines the behaviour of 
the masses. We expect this mechanism to apply to high- 
spin mesons as well. Therefore we expect that for a me- 
son of fixed spin and low excitation number, the mass 
will decrease as a function of temperature. We will in- 
deed find this behaviour. However, there is now an ad- 
ditional interesting phenomenon, which only appears for 
high-spin mesons, since their spin can vary significantly. 
Spin is a parameter which controls the length of the string 
as well the distance of the tip of the string to the hori- 
zon. Hence, a natural question arises as to whether the 
confinement-deconfinement transition of mesons (i.e. the 
process of "dissociation" of the mesons), takes place at 
the same temperature for all mesons, irrespective of their 
masses and spins? Intuitively, one expects that higher- 
spin mesons will dissociate first, given that they "hang" 
closer to the horizon. The question is then what is the 
critical value of the spin as a function of temperature? 
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0.25 0.5 0.75 1 1.25 1.5 1.75 2 

FIG. 6: The boundary curve (|24[) . Rotating strings have to lie 
above this curve in order for their action to be real. The var- 
ious curves correspond to various values of the frequency ui. 
The horizon is located at ut = 2. These curves are indepen- 
dent of the probe brane location. 

To answer the questions posed in the previous para- 
graph, let us consider large spinning mesons at finite 
temperature. These are modelled by rotating U-shape 
strings hanging from the probe brane in the middle pic- 
ture of figure [TJ The relevant part of the background 
metric @ is 




We go to the static gauge for the string action and make 
the following ansatz for the rotating configuration, 4 

t = t , p = a , u = u(p) , ip = ujt . (22) 

We see that the ansatz has the the same form as in the 
zero-temperature case [3l|. Hence, just like for the Wil- 
son loop, the only effect of finite temperature will be that 
the details of the shape u{a) change as the temperature 
is increased. 

One could wonder why such an unmodified ansatz is 
correct, given that the string is moving with angular ve- 
locity, and velocity dependence has been shown to mod- 
ify the zero-temperature ansatz for vertical strings due to 
"dragging" effects [0, QJ, OS 0, US ES El • To answer 
this question we need to see if the ansatz ([22]) allows for 
physical solutions. If it does, this would imply that an- 
gular velocity of a meson in a hot medium does not lead 
to the distortion of rotation, i.e. one does not need to 



4 When solving the equations of motion numerically, the fact that 
this gauge leads to a dti/dcr-derivative which blows up at the 
flavour brane is rather inconvenient. For numerical purposes, a 
better gauge choice is, for instance, given by p + u = a. 
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FIG. 7: Comparison of numerical solutions (solid curves) with the boundary curve (|24l) (dashes), for lu = 0.5 and u = 1 
respectively. Although the numerical solution lies manifestly above the boundary curve, as it should, it deviates quite a bit 
from it. The numerical solutions where obtained by fixing u(0) = ut + e for a small e, as well as w'(0) = 0. 



apply a force to maintain the rigid rotation: an initially 
rigidly rotating meson continues to rotate rigidly. This of 
course does not mean that the u(p) profile is unaffected 
by the temperature. A separate question is what is the 
effect of linear motion of this colour singlet. This will be 
explored in section HVBI 

To answer the first question, about the possibility of 
having a rigidly rotating meson in a quark gluon plasma, 
let us compute the action for the meson (122p . This ansatz 
leads to 



/ 



S= drdp 




W 2 i%4 

f(u) u 3 



(/(«) - P 2 UJ 2 ) 

(23) 

Before we analyse the resulting equations of motion, let 
us first make a few general remarks. We see that posi- 
tivity of the argument of the square root in (|23p requires 
that f(u) > p 2 ui 2 . This in turn means that for a given 
angular frequency lj, the string solution u(p) has to lie 
above the curve 



u{p) > 



(1 



^2)1/3 



(24) 



In figure [S] we have depicted these curves for various val- 
ues of lj. Firstly, we see that the tip of any string (the 
p = point) is in principle allowed to touch the horizon 
for any angular frequency u. Secondly, for given u>, the 
maximal allowed extent of the string is determined by 
the intersection of the curve with uq, and is given by 



_ 1 I U T 
Pmax — \ / -I 

Uq 



(25) 



The equations of motion following from the action 
are given by 



I u' 



(f(u)-pW) 



•/'(«) 



/2 2 2 



/(«) 



xl D4 



if 



'\2 



GoT-(p')W =0, (26) 



3u 2 
R 3 



where ^/TTT is the density of Nambu-Goto action (|23|) . 
One can verify that varying the action first and then in- 
serting the ansatz yields the same result. In addition 
we have to take care of the boundary terms. If one 
were not yet in the static gauge, one would impose Neu- 
mann boundary condition for the p direction and Dirich- 
let in the it-direction. These two, when combined imply 
that string end "orthogonally" on the brane worldvolume 
du/dp — > oo. 

Finally, the expressions for the energy and the angular 
momentum carried by the string are given by 



E 



dcr. 



1 



( u 



\Rji4, 



f(u)(p') 2 + u'- 



J= da 



1 



LUp 



(27) 



(28) 



where 
action 



RmJ f(u) 
is once more the density of the Nambu-Goto 



(One could have hoped that the boundary curve 
becomes a good approximation to the real string solution 
in some regime. However, this does not seem to be the 
case, as is easily seen from a comparison of this curve 
with some real solutions, figure [7]) 



A. Properties of the spectrum and melting of 
large-spin meson 



Due to the complexity of equation (|26p . we will anal- 
yse the spectrum of large spinning strings numerically. 
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FIG. 8: Shapes of rotating strings with fixed position of the 
flavour brane «o = 20 (i.e. almost fixed L/R, see figure [2} 
and fixed temperature Ut = 1, for various values of the fre- 
quency LU. 



Nevertheless, interesting qualitative features can easily 
be extracted from this analysis. 

Let us first discuss the spectrum for fixed temperature. 
By analysing the shapes of the string for various values 
of to, we see that as lu decreases, i.e. the spin of the meson 
increases, the string endpoints get more and more sepa- 
rated, and the U-shaped string penetrates deeper to the 
horizon. It becomes more and more rectangular shaped. 
A few characteristic shapes, for mesons with fixed T and 
m q but varying lu (or equivalently J) are depicted in fig- 
ure El 

Integration of the equation yields the dispersion 
relation in parametric form, i.e E(u>) and J(lu) at fixed 
temperature. These curves, together with similar curves 
at zero temperature, are shown in the plots in figure [9] 
We see that at zero temperature, the spin increases as the 
angular frequency lu is decreased (just like in flat space, 
where this behaviour is essentially determined by the con- 
dition luL = 1 which says that the endpoints move with 
the velocity of light). As the spin increases, the string 
hangs deeper into the geometry and hence the energy of 
the meson goes up. The same happens at non-zero tem- 
perature for large values of lu, i.e. when the string is short 
and thus far away from the horizon. However, as lu is de- 
creased, even though the length of the string continues 
to grow, the effective tension of the parts of string in the 
region near the horizon decreases rapidly. This effects 
leads to a slower growth of spin and energy, and hence 
to the appearance of a maximum in the energy and spin. 
This is shown in the solid (red) curves of the plots in 
figure [51 

When we eliminate the variable cj from the two plots in 
figure^ to obtain a plot of the energy versus the spin, we 
get figure [10] The appearance of a maximum in energy 
and spin is again clearly visible. We also see that the two 
states with identical spin J actually have different energy: 
the ones with smaller lu are more energetic than the ones 



with larger lu. Therefore, the upper branch to the left of 
the maxima in the plots of figure[5]is presumably unstable 
and will decay to the lower branch. 

Apart from the two configurations in figure [TUI with E\ 
and E2 for fixed J, there is yet another configuration one 
may consider; this configuration is in spirit closer to the 
dissociation of Wilson loops. Namely, one could consider 
the configuration of two dragged single quarks (to be de- 
scribed in more detail in the next section), with linear 
velocities v which are such that the whole system carries 
the same angular momentum as our mesons. Then, one 
could compare the energy of such a configuration with 
the energy of the meson before dissociation. However, 
as shown in this configuration always carries infinite 
energy, and thus cannot be the one which acts as initial 
data for the dissociated meson. 5 

The plot of J versus lu = w(ubottom) is similar to the 
plot of the inter-quark distance versus the minimum dis- 
tance to the horizon found for Wilson loops. For Wil- 
son loops, it was found that there is a maximum inter- 
quark distance. For spinning strings, such a maximum 
distance appears too. However, this maximum distance 
is not reached for the same value of lu as the value for 
which the maximum spin is attained. A priori, this is a 
source of worry, since it could mean that the relation be- 
tween J and the inter-quark distance L is not monotonic. 
However, we have verified numerically that the maximum 
inter-quark distance always occurs in the unphysical part 
of the spectrum; the relation between J and L is mono- 
tonic in the physical regime (this is somewhat similar to 
the monotonicity of L versus yx depicted in figure ^ . 

From the analysis above, we see that for a fixed tem- 
perature there is a maximal value of the spin which a 
meson can carry. It is natural to interpret the tempera- 
ture at which this happens as the critical temperature at 
which a meson of spin J max melts. Thus we conclude that 
(as intuitively expected) the dissociation temperature of 
large-spin mesons is spin dependent. As the temperature 
increases, the maximal value of the spin that a meson can 
carry decreases (see figure [TTj) . i.e. for given quark mass, 
higher-spin mesons melt at lower temperature. We also 
see that for meson of fixed angular momentum, the en- 
ergy decreases as a function of temperature, see figure 
This is the same behaviour as was observed for the low- 
spin mesons. However, the four-dimensional size of the 
meson (i.e. the distance between the two endpoints of 
the U-shaped string measured along the probe-brane) in- 
creases as the temperature goes up; see figure [T51 Finally, 
by considering the shapes of strings at different temper- 
atures (figure Q3]) we observe that, as the temperature is 
increased, the largest allowed strings become less and less 
rectangular shaped, and more and more U-shaped, thus 
leading to the breakdown of the four-dimensional picture 



5 In this sense the dissociation of spinning strings is different from 
the dissociation of Wilson loops [32] , |33[ . 
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FIG. 9: The energy (squared) and angular momentum of U-shaped strings, as a function of the frequency ui. The dashed (blue) 
curves are for zero temperature, while the solid (red) curves are for the intermediate-temperature regime. 




FIG. 10: Left panel: Energy versus angular momentum for L/R = 0.24 in the intermediate-temperature phase, T/T c = V%- 
The lower segment (solid) is for w > 1, the upper one (dashed) for uo < 1. There is a maximum reachable spin. Right panel: 
Similar plot, but now for two different temperatures. The blue (upper) curves are for a lower temperature T/T c — 1, the red 
(lower) curves for a higher temperature T/T c = \/2. Observe that for a fixed spin, the energies go down as the temperature 
goes up. 



of a string with two massive, point-like strings attached 
to it [3lj . This is an important observation if one wants 
to construct four-dimensional effective string models for 
mesons at finite temperature. 



B. Drag effects for mesons 

It was recently observed that the drag force experi- 
enced by free quarks moving through a hotgluon plasma 
can be modelled in string theory [ll|, [l2|, [ljl, 0, EH, [H, 
[TtJ ■ In this setup, the quark is a single string ending with 
one end on the flavour brane and extending all the way 
down to the horizon. Because a strictly vertical string 
moving rigidly through the background would not have 
a real action (j2"5)) . the string has to be bent when it is 
"pushed" through the plasma. In addition the bent string 
does not end anymore orthogonally on the brane. This 



means that one has to apply a force on the string end- 
point, or in other words, one has to "drag" the string in 
order to keep it moving. 

In more detail, this conclusion follows from considering 
the string action evaluated on the "naive" ansatz describ- 
ing a rigid straight string moving with velocity v x , 

t = T, u = a , x = v x t . (29) 

Plugging this ansatz into the action (f2"5| , we get 



Squark ~ \T /S)' (30) 

For this expression to be real, a purely vertical string 
cannot stretch lower than 

u>u ^ "° = fl ^i/s > ( 31 ) 
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FIG. 11: Analysis of maximum spin for U-shaped strings 
in the intermediate-temperature regime, as a function of 
temperature, for fixed L/R = 0.24. Angular momentum 
is normalised as in (|28[1 . There is of course a divergence 
as T/T c — > 0, but before that point is reached, one encoun- 
ters the phase transition to the low-temperature background 
(vertical dotted line at T/T c = 1). 
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FIG. 12: Analysis of the temperature dependence of the en- 
ergy of U-shaped strings. Displayed are the squared energies 
of a state of spin J ~ 1.5, which dissociates at T/T c — 3.16 
when L/R = 0.24. The diagonal dotted curve shows the 
maximum possible energy for a given temperature. Note the 
qualitative similarity with the low-spin plot, figure [5] 



and it clearly cannot stop there when there is no sec- 
ond flavour brane present. The true shape is obtained 
by adding a c-dependent term to the ansatz for x. The 
string now bends as it reaches down to the horizon. The 
bent, moving string does not, however, satisfy Neumann 
nor Dirichlet boundary conditions, neither at the horizon 
nor at the probe brane. Hence the world-sheet momen- 
tum does not vanish at the string endpoints: there is a 
world-sheet momentum flow into the horizon, which has 
to be "added" by a force dragging the string on the up- 
per, probe brane side. This effect translates to a drag 
force in the gauge theory picture. 

Hence in summary, we see that there are two effects 
happening as one tries to move a single string in the 
hot background: firstly, the string shape is modified in 



FIG. 13: The four-dimensional size of a large-spin meson of 
fixed spin as a function of T/T c : when the temperature goes 
up, the meson expands. The dotted curve denotes the size of 
the maximum-spin meson at a given temperature. 



a way which depends on the temperature and velocity, 
and secondly, in order to preserve the motion one needs 
to apply a force. 

It is a natural question to ask if such effects are also 
exerted on spinning mesons, as the temperature is turned 
on. One might expect that even a simple rotating motion 
experiences a drag effect, but the analysis of the previ- 
ous subsection shows that this effect docs not exist. In- 
stead, the rotating string is always sufficiently high above 
the curve beyond which the action would turn imaginary. 
No world-sheet momentum gets lost behind the horizon. 
On the other hand, the first effect is still present. The 
bending of the rotating meson does depend both on the 
angular velocity and on the temperature. 

However, there is another type of motion which one can 
consider, and which is more similar to the moving-quark 
situation. Namely, one can consider a linear motion of 
the meson in a direction orthogonal to the plane of rota- 
tion. 6 A suitable ansatz for this motion is 



t = T , p = a , U = U{p) , ip = LOT , y = v y T . 



In this case the action becomes 



(32) 



5' 



dr dp * 



R 3 



/(«) 



(/(«) 



p 2 L0 2 



(33) 

The only modification with respect to the rotating meson 



is the addition of a term 



to the last factor under 



the square root. The condition for the action to be real 



6 One could also consider motion in the direction in which the 
string rotates. In this case the equations of motion are much 
more involved. We do, however, expect that the conclusion 
reached by studying this motion will be similar to those reached 
by studying motion in a transverse direction. 
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FIG. 15: Analysis of the effect of a transverse velocity on spinning U-shaped strings. The left panel depicts the relation 
between J and u> for fixed position of the flavour brane. The upper curve is for v y = 0, while the other curves display results for 
increasing values v y — 0.2, 0.4, 0.6, 0.8, 0.9, 0.98. The panel on the right shows the effect of a transverse velocity on the shape, 
keeping the quark masses and spin fixed (corresponding to the dashed black line in the left panel). 



u > 



(1 



p 2 w 2 



V 2)l/3 



(34) 



Hence, if the ansatz (f3"2"j) makes sense, there should exist a 
solution to the equations of motion such that the spinning 
configuration lies entirely above this curve. If that is not 
the case, the ansatz has to be modified. 

We have verified that the equations of motion obtained 
from (|33|) are the same as those of the full system after 
insertion of the ansatz (|32[) . Subsequently, we have in- 
tegrated the equations of motion using the same tech- 
niques as for the zero-velocity case. The results are de- 
picted in figure [TSl The left panel shows that, as the 
transverse speed is increased, the maximum possible spin 
decreases. If we take a state of fixed spin J and look 
at how its shape changes as we increase the transverse 
speed, we obtain the second panel in figure [T51 (this cor- 
responds to following the states which lie on the black 



dotted curve in the first panel). All these configurations, 
up to the maximum-spin ones, lie safely above the limit- 
ing curve (|3~4"]) . Note how the size of a fixed-spin meson 
increases slightly as the transverse velocity is turned on. 

Thus, the important physical conclusion seems to be 
that colour singlet states do not experience any drag ef- 
fect. From the field-theoretic point of view this makes 
sense, as the absence of a colour monopole moment means 
that these states do not couple directly to gluon degrees 
of freedom. Because of the absence of drag, mesons do 
not experience any energy loss when propagating through 
the quark-gluon plasma: no force is necessary to keep 
them moving with fixed velocity. 7 Note moreover that 



7 We should emphasise that what we have computed here is dif- 
ferent from the computation of |17j |. In that paper, a situa- 
tion is considered in which two quarks, being endpoints of a 
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FIG. 16: The behaviour of the four-dimensional size L for 
maximum-spin mesons as a function of the transverse velocity 
(for a particular choice of m q and T = T c ). The blue curve 
depicts L 4d (v = 0) ■ (1 — v 2 ) 1 ^ 4 , as obtained analytically for 
Wilson loops in [341 ]. 



the moving meson does not need to bend in the direction 
of motion, i.e. there is no y(u) term present in (|32p . How- 
ever, the shape of the meson in the (p, u) plane certainly 
is velocity dependent. As shown on the right panel of 
figure IT5l the distance to the horizon decreases, leading 
to a lower melting temperature. 

Figure [T5] (left panel) also shows that the maximum 
spin goes down as the velocity goes up. There is thus a 
critical velocity beyond which a meson of fixed spin has 
to dissociate. Similarly, the four-dimensional size of the 
meson decreases with increasing velocity, as can be seen 
in figure Uni The data is approximated rather well by the 
relation 

-^max-spin ( v ) ~ ^max-spin ( v = 0) • ( 1 ™ l> 2 ) ^ . (35) 

This fit was motivated by the analytic results of [341 ] . in 
which a similar dependence on the velocity was found 
for the screening length, or more precisely, the maximum 
interquark distance for Wilson loops in an AdS black hole 
background. See also related results in [35l |. 

Finally, one should keep in mind that our analysis in 
this section is purely classical, and restricted to the sigma 
model, without any coupling to the fields on brane. In 
principle there will be classical (and quantum) energy 
loss due to radiation effects on brane, because the meson 
is essentially a dipole. However, the interpretation of this 
energy loss is different, as it corresponds to the decay of 
high-spin mesons into low-spin ones [36l . 1371 . |38| , both at 
zero and at non-zero temperature. 



"twisted" Wilson loop, move away from each other. The re- 
quired energy influx is higher at finite temperature, giving rise 
to a four-dimensional drag force. For us, the two quarks which 
make up the meson move in the same direction, and we do not 
change the U-shaped string as a function of time. 



V. SPECTRUM AT HIGH TEMPERATURE 

Finally, let us briefly discuss the high-temperature 
phase. In this phase, the minimal free energy of the sys- 
tem is attained when the two stacks of branes are dis- 
connected (see the third panel of figure [J). The profile 
of the left and right stacks of branes is characterised by 
u' = dw/dx 4 — > oo and the induced metric on the probe 
branes and probe anti-branes takes the form 




i.e. it is the same as the background metric without the 
544 (dx 4 ) 2 term. Considering now the action of the gauge 
fields on each of the two stacks of branes, it is straight- 
forward to see that the action will be similar to the one 
at intermediate temperature, but with 7 = / . The 
differential equation for the modes is now 

- u 1 ' 2 /(«) d u (u 5 / 2 f(u)d u ^ (n) ) = Rl 4 m 2 n V> (n) . (37) 

The normalisation conditions are similarly simplified, 




The first thing to observe now is that the massless pion 
disappears from the spectrum. This mode, which would 
be given by 4>^ — u _5 ' 2 /(it) _1 , is no longer normalis- 
ablc. Computation of its norm leads to the integral 



f* 00 

/ duu 5/2 









which, while convergent at the upper boundary, is di- 
vergent at the lower boundary because f(u) ~ \Ju — ut 
for u ~ ut- Here, it is crucial that the D8-branes ex- 
tend all the way down to the horizon at u — ut (recall 
that at finite temperature, these /-factors did not lead 
to a divergence because the integral had a lower limit 
at u = uq). In accordance with the fact that chiral sym- 
metry is restored in the high-temperature phase, we see 
that the Goldstone boson has disappeared. 

In order to see whether the remainder of the spectrum 
is discrete or continuous, we write our equation (|37|) as 
(we set i?D4 = ut = 1 from now) 

f(u) d u (u 5 ' 2 f{u)d u il> {n) ) + u- 1 ' 2 m 2 W) = . (40) 

If we expand this near u = 1 we get 

9(u-l)dJ(u-l)d u ip(ufj +m 2 ijj(u) = 0. (41) 
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The general solution to this equation is 

— log(tt - 1)J + c 2 sin y— \og(u - 1)J . 

(42) 

These modes are normalisabile using the norm (f38|) for 
any value of the constants c\ and c^- 

If one now shoots from infinity, starting with the nor- 
malisable mode for a fixed m 2 , one can always match 
this to a regular solution at the horizon because there 
are two free parameters available there. Because there 
are two regular solutions at the horizon the spectrum 
will be continuous in the high-temperature phase. 

Note that this type of analysis could in principle be 
made also for the zero- and intermediate-temperature 
phases. In those cases, one will find that again both solu- 
tions are regular at the tip of the probe brane. However, 
there is now an additional condition, namely that the so- 
lution can be continued to the other side of the probe 
brane and remain normalisable. So in these two phases, 
it is not regularity at the tip, but rather the condition 
that the function is either odd or even which makes the 
spectrum discrete. 



Note added 

After publication of the first version of this paper, (34| 
appeared, which contains a computation of the screening 
length at finite temperature (the maximum interquark 
distance for Wilson loops). Our results leading to fig- 
ure [15] agree with the analytic behaviour ^max-spinC 1 ") ~ 
^max-s P in( v = 0) (1 - u 2 ) 1 / 4 found there; we have added 
figure [16] and the paragraph around (|35|) to make this 
more explicit. 
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